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We review recent developments in the QCD description of jet production and modifi- 
cation in reactions with heavy nuclei at relativistic energies. Our goal is to formulate a 
perturbative expansion in the presence of nuclear matter that allows to systematically im- 
prove the accuracy of the theoretical predictions. As an example, we present calculations 
of inclusive jet cross sections at RHIC, Z°/7*-tagged jet cross sections at the LHC, and jet 
shapes that include both next-to-leading order perturbative effects and the effects of the 
nuclear medium. 

§1. Introduction 

In the past several years important developments in jet finding algorithms 11 com- 
bined with advances in detector technology and experimental analysis have enabled, 
for the first time, measurements of jets® in nuclear reactions at very high energies 
at the Relativistic Heavy Ion Collider (RHIC).® Observation of jets is among the 
first physics results reported at the Large Hadron Collider (LHC), 4 ' which also has 
an active heavy-ion program.® In nucleus-nucleus (A+A) collisions, jet observables 
are expected to be much more discriminating® with respect to the underlying modes 
of parton propagation and energy loss in strongly-interacting matter® than existing 
studies of leading particles and leading particle correlations. Developments if theory 
are on the way to guide the heavy- ion experimental jet programs at RHIC and at 
the LHC and to help interpret current and upcoming physics results. 

To take full advantage of jet physics, calculations at next-to-leading order (NLO) 
and beyond in perturbative Quantum Chromodynamics (QCD) are required.®'® In 
heavy-ion reactions there is the added complication of a soft background medium 
that these jets must traverse - the quark-gluon plasma (QGP). There are multiple 
ways to describe its properties that include but are not limited to temperature T, 
coupling strength g med ', density p, energy density e, parton rapidity density dN 9 /dy, 
Debye screening scale mo, parton mean free paths X q , X g , transport coefficient q, 
and the strength of the background gluon field (F + - L F^). Relations between some 
of these quantities can be derived for specific models of the nuclear matter. In 
practice, however, one always needs independent constrains from experimental data 
or lattice simulations since the QGP that A+A reactions aim to produce may be 
strongly-coupled or non-perturbative. 

The production of jets and the medium-induced bremsstrahlung at scales Q 2 ~ 
E\ Aq CD and Q 2 ~ (gT)Er S> Aq CD , on the other hand, can be treated pertur- 
batively. In fact, parton energy loss processes in the QGP and one-loop perturbative 
corrections are formally manifested in experimental observables at the same order 
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Table I. Perturbative expansion relevant to large Q 2 processes in the presence of a nuclear medium. 



0(atg), 0(a 2 af eA -). It should be noted, however, that the values of a s and a™ ed ' 
can differ considerably. The emerging perturbative expansion for jet production and 
modification in heavy-ion reactions can easily be generalized to higher orders and 
is shown in Table [H Self-consistent fixed-order theoretical calculations are the ones 
that take into consideration all terms along the counter diagonals of this matrix. At 
present, only inclusive a™ ed ' induced gluon bremsstrahlung has been evaluated and 
the corresponding next-to-leading order jet calculations are presented below. 



§2. Inclusive jet production at RHIC 



At NLO, the inclusive and tagged jets cross sections can be expressed schemat- 
ically as follows: 



d<T je t(+t ag ) = \ { da[2 -> 2]S 2 ({p,yAh) + j { do[2 -> 3]5 3 ({p, y, 0} 3 ) • (2-1) 



Here, pi, yi, fa are the transverse momentum, rapidity, and azimuthal angle of the i-th 
particle (i = 1,2,3), respectively, and a[2 — > 2], a[2 — > 3] represent the production 
cross sections with two and three final-state partons. S2 and S3 are phase space 
constraints and S2 = Yli=i = Yld=i <5(-Et; — ET)${yi — y) identifies the jet with its 
parent parton. Hence, it is only at next-to-leading order that the dependence of the 
experimental observables on the jet cone radius R, the jet finding algorithm, or the 
trigger particle energy can be theoretically investigated. For an angular separation 
Rij = sj {yi — yj) 2 + (fa — <pj) 2 , defined for any possible parton pair (i, j), 

S 3 = ]T 8{E Ti + E T] - E T )5 ( ET E Vl + E E m - y] 9 (R tJ < R rc ) (2-2) 

i<j 1 i 

+ e s ® n 9 > {ET )t ET ^ ) > r ™ = min ( r «* r > ET ^ Er s s r ) ■ 

i J-ti V max{E Ti ,E Tj )J \ max{E Tt ,E Tj ) J 

In Eq. ()2-2p R IC determines when two partons should be recombined in a jet. Here, 
1 < Rsep < 2 is introduced to take into account features of experimental cone algo- 
rithms, employed to improve infrared safety. Eq. (|2-2p establishes a correspondence 
between the commonly used jet finders and the perturbative calculations to 0(a^ s ) 
with the goal of providing accurate predictions for comparison to data. For example, 
R sep = 2 yields a midpoint cone algorithm and R se p = 1 corresponds to the kx algo- 
rithm.® 'E2> R is the cone size or parton separation parameter, respectively. We have 
compared NLO calculation® of the inclusive jet cross section at ^fs = 200 GeV p+p 
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Fig. 1. Comparison of the NLO calculation to STAR experimental data^ on the inclusive jet cross 
section for R=0.4 (left panel). The variation of the jet cross section with the cone size R for 
Et = 10, 40 GeV around midrapidity at RHIC is also shown (right panel). 



collisions at RHIC to the STAR experimental measurement which uses a midpoint 
cone algorithm 11 ' in the pseudorapidity range 0.2 < rj < 0.8. Very good agreement 
between data and theory is achieved with a standard choice for the renormalization 
and factorization scales fJ>R = (J-f = EtW^ as shown in the left panel of Figure [TJ 
Variation of these scales within (Et/2,2Et) leads to less than (+10%, —20%) varia- 
tion of the jet cross section. We also found that there is a significant dependence of 
da }et I dydET on the cone size R, which, even in p+p reactions, can exceed a factor 
of two. This is illustrated in the right panel of Figure [TJ Analytically, the ln(i?/i?o) 
scaling of the cross section can be understood from the 1/r angular behavior of the 
perturbative QCD splitting kernel at high energies. 

When compared to a parton shower in the vacuum, the medium-induced quark 
and gluon splittings have noticeably different angular and lightcone momentum frac- 
tion dependencies P>E3 In particular, for energetic partons propagating in hot and 
dense QCD matter, the origin of the coherent suppression of their radiative energy 
loss, known as the Landau-Pomeranchuk-Migdal effect, can be traced to the cancel- 
lation of the collinear radiation at r < mc/(w(mfl, X g ,E T )). 13} Here, the Debye 
screening scale ~ gT and (cj) ~ few GeV. Thus, the medium-induced com- 
ponent of the jet, which is given by the properly normalized gluon bremsstrahlung 
intensity spectrum ^ med (r, R) tx dl ra,d /dujdr within the cone, has a characteristic 
large-angle distribution away from the jet axis. This is illustrated in Figure [2] for 
central Au+Au and central Cu+Cu collisions at RHIC. We emphasize that accurate 
numerical simulations, taking into account the geometry of the heavy-ion reaction, 
the longitudinal Bjorken expansion of the QGP, and the constraints imposed by its 
experimentally measured entropy density per unit rapidity,^ 1 have been performed 
for all physics results quoted here. 

One can exploit the differences between the vacuum and the in-medium parton 
showers by varying the cone radius R (Rijet < R) and a cut p™ in (E^ > p™ in ) for 
the particles "i" that constitute the jet, to gain sensitivity to the properties of the 
QGP and of the mechanisms of parton energy loss in hot and dense QCD matterP 
This is illustrated in the insert of Figure [2j The most easily accessible experimental 
feature of jet production in nuclear collisions is, arguably, the suppression of the 
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Fig. 2. The differential jet shape in vacuum ifj vac ' (r, R) is contrasted to the medium-induced 
contribution i/) mcd '(r, R) by a Et = 30 GeV quark in Au+Au and Cu+Cu collisions at 
y/SNN = 200 GeV. The insert illustrates a method for studying the characteristics of these 
parton showers (left panel). Transverse energy dependent nuclear modification factor R 3 ^ A f° r 
different cone radii R in b — 3 fm Au+Au collisions. Inserts show ratios of jet cross sections for 
different R in nuclear reactions versus Et (right panel). 



inclusive cross section in heavy-ion reaction compared to the binary collision scaled, 
oc (-/Vbm), production rate in elementary nucleon-nucleon reactions!^* 

let/ mim da AA (E T ; R, p^ in ) I da pp (E T ; R, p™ in ) . . 

Rf A (Er;R,Pr) = j (AW . (2.3) 

Eq. (|2-3p defines a two dimensional jet attenuation pattern versus R and p™ in for ev- 
ery fixed Et- In contrast, for the same Et, inclusive particle quenching is represented 
by a single value related to the R — > and p™ in S> ( uj ) limit in Eq. (|2-3p . Thus, jet 
observables are much more differential and, hence, immensely more powerful than 
leading particles and leading particle correlations in their ability to discriminate be- 
tween the competing physics mechanisms of quark and gluon energy loss in dense 
QCD matter and between theoretical model approximations to parton dynamics in 
the QGP. 

We calculate the medium- modified jet cross section per binary nucleon-nucleon 
scattering as follows (p^ m = 0): 

1 da AA (R) _ yi ^ 1 da?r^°(R) (24) 

(N hm ) dyd?E T ~ J e=0 ^' 9(e ' ^ (1 - (1 - f q , g ) ■ e) 2 dyd*E> T ' { ^ } 

Here, -P g , g (e, E) is the probability distribution for the parent quarks and gluons to lose 
a fraction e = Y^i Wi/E of their energy due to multiple gluon emission in the QGP. In 
Eq. ([231) d<j^j> Nh °{R)/dyd 2 E T is the differential cross section which includes the 
known cold nuclear matter effects^ and (1 — f q>g ) ■ e represents the fraction of the en- 
ergy of the parent parton that the medium re-distributes outside of the cone of radius 
R. The measured cross section is then a probabilistic superposition of the cross sec- 
tions of protojets of initially larger energy E' T = £t/(1 — (1 — Our results for 
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the nuclear modification factor of inclusive jets R?2 A in central Au+Au collisions with 
■s/snn = 200 GeV at RHIC are presented in the right panel of Figure [2j Experimen- 
tal data on leading ir° suppression for these reactions is only included for reference. 
A continuous variation of R^aa with the cone radius R is clearly observed and shows 
the sensitivity of the inclusive jet cross section in high-energy nuclear collisions to 
the characteristics of QGP-induced parton shower. For R < 0.2 the quenching of jets 
approximates the already observed suppression in the production rate of inclusive 
high-pT particles. It should be noted that in our theoretical calculation CNM effects 
contribute close to 1/2 of the observed attenuation for Ej< > 30 GeV. These can be 
dramatically reduced at all Et by taking the ratio of two differential cross section 
measurements for different cone radii R\ and R 2 , as shown in the insert. 

§3. Z°/7*-tagged jets at the LHC 

We begin by discussing the cross section for Z°/7*-tagged jet production in 
p+p collisions. It is instructive to first consider the leading order (LO) result, from 
which one can understand the underlying production processes and appreciate why 
the Z° boson was originally considered as a suitable tag for the initial associated jet 
energy PS In the collinear factorization approach, the Z°/7*+jet cross section reads: 

d°_ _ f(xi,n)f(x 2 , ji)\M\ 2 _ ^ 

dy(z)dyQ et )d 2 p T{z) d 2 p T Q ct) ^ (2vr) 2 4xi x 2 S 2 T{z) T(jet) 

Here, f(xi, /x) are the parton distribution functions and \M\ 2 are the relevant squared 
matrix elements. The 5-function constraint on the transverse momentum of the jet 
is valid only at this order and only at the partonic level. 

The main advantage of the next-to-leading order Z°/7*+jet+X calculation® is 
the ability to precisely predict the transverse momentum distribution of jets associ- 
ated with a dimuon tag in a narrow p? interval.^ 1 Beyond tree level, the momentum 
constraint that the Z° boson measurement provides is compromised by parton split- 
ting and Z-strahlung processes. We demonstrate this in the left panel of Figure O 
which shows the differential cross section for jets tagged with Z°/j* — > fi + pT in p+p 
collisions at ^/snn = 4 TeV. We implement acceptance cuts of \y\ < 2.5 for both jets 
and final-state muons, and, as mentioned above, constrain the invariant mass of the 
dimuon pair to the interval M z ± 3r z , where M z = 91.1876 GeV and F z = 2.4952 
GeV. This is the kinematic acceptance range in which we will evaluate all results 
that follow and which can be easily adjusted to match upcoming experimental mea- 
surements. For the cross section shown in Figure [3j the tagging Z°/7* is required 
to have 92.5 GeV < pt < 112.5 GeV. The LO result restricts the pr of the jet to lie 
exactly within this interval, consistent with Eq. (|3Tp . As seen in the left panel of 
Figure[3l at NLO the deviations from this naive relation are very significant. We have 
included results for three different values of the jet cone radius, R = 0.2, 0.4, 0.8. 
The variation of the cross section with R around PTQet) ~ Pt(Z) arises from the 
interplay between the amount of energy that is contained in the jet and the number 
of reconstructed jets (one or two). For Pr(jet) ^> Pt(z) or PTQet) ^ Pt(z) the two 
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Fig. 3. Transverse momentum distributions of a jet associated with Z°/ r y*(^ /i + + ^T) tag to 
0(Gfois) and 0(Gfo?s) for transverse momentum cut 92.5 GeV < pt < 112.5 GeV on the 
tagging particle (left panel). The NLO pT-differential cross section per nucleon pair for these 
tagged jets in central Pb+Pb collisions at the LHC is also presented for a cone radius R = 0.2 
(right panel). The ratio of the two cross sections is shown in the insert. 



final-state partons are well-separated and identified as different jets. The falloff of 
the differential cross section relative to its peak value at pj> Q c t) = Pt (z) is then con- 
trolled by the QCD splitting kernel (the part related to the large lightcone parton 
momentum) and there is no dependence on the cone radius. 

In order to quantify the inability of the Z /^* tag to constrain the momentum 
of the jet we calculate the mean p? = {pt Qct)) an d standard deviation Ap T Q et \ = 

y^TQet)) ~~ ^T(jet)) 2 f° r each of the curves in Figure [3l Our results are presented 
in Table 2. The standard deviation for the LO curves is not strictly zero because 
of the finite pr width of the tagging Z°/j* interval. The NLO curves exhibit a 
similar (px(jet)) as ^ ne ^0 result, with (pTQet)} increasing as the cone radius in- 
creases. However, in going from LO to NLO there is a significant jump in Ap^^y 
The width of the jet momentum distribution quadruples for an energetic tag. The 
very large values of Ap T Q e t) / (Pt (jet)) ~ 25% at NLO create serious complications 
for experimentally tagging the initial associated jet energy in both p+p and A+A 
collisions. While additional cuts can be considered, such as the requirement for a 
single jet within the experimental acceptance that is exactly opposite the tagging 
particle in azimuth, these will reduce the already small projected multiplicity for this 
final state in heavy-ion reactions. 

Accounting for the fact that the Z°/j* and their decay dileptons escape the 
region of dense nuclear matter unaffected by the strong interaction, the modified jet 
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Table II. Mean pr and A{pr) for Z°/7*-tagged jets at the LHC from LO and NLO calculation. 



The theory of jets in nuclear collisions 



7 



cross section can be calculated as follows: 

= vf rff *jM daq ' 9 ( p Q 

d 2 V {Z) cPv Q [1 - (1 - f q , g )e] 2 d 2 P ( z)d 2 P iict) \[1-(1- f M )e)] 

Eq. (|3-2p resembles the result for inclusive jet suppression and we have integrated over 
yjet and in finite rapidity intervals. Our results are presented in the right panel of 
Figure El where the tagging Z°/~/* is required to have 92.5 GeV < px < 112.5 GeV. 
The main physics effect that this figure illustrates is the QGP-induced modification 
to the vacuum parton shower. Specifically, its broadening^ 1 implies that part of the 
jet energy is redistributed outside of the jet cone and the differential jet distribution 
is downshifted toward smaller transverse momenta. The smaller the jet cone radius 
the more pronounced this effect is. As the jet cone radius is increased, the medium- 
modified curves approach the p+p result, as more and more of the medium-induced 
bremsstrahlung is recovered in the jet. 

A generalization of the cross section ratios per binary nucleon-nucleon collisions 
for tagged jetsf^ 1 Ij^a, is given in the insert in Figure [3l The variation in the 
magnitude of Pft, at transverse momenta larger than the px of the tag is controlled 
by the shape of the jet spectrum. The most striking feature is the sharp transition 
from tagged jet suppression above px (z) to tagged jet enhancement below px <z) ■ For 
the example shown in the left panel of Figure O the enhancement can be as large as a 
factor of ten. This transition from suppression to enhancement is a unique prediction 
of jet quenching for tagged jets and will constitute an unambiguous experimental 
evidence for strong final-state interactions and parton energy loss in the QGP. 



§4. Jet shapes at RHIC and at the LHC 



Integral and differential jet shapes, defined as: 

were historically the first observables used to study jet sub-structure and intra-jet 
energy flowP To <D{al) 4>(r,R) can be evaluated from the LO parton splitting 
functions in perturbative QCD^" and the analytic results can be generalized to allow 
for a minimum particle or calorimeter tower energy p™ m in the definition of the jetP 1 
Such experimental cuts can be particularly useful in reducing the large background 
of soft particles in the high multiplicity environment of heavy-ion collisions. 
Analytically, jet shapes are evaluated as follows^ ® 

V, vac - (r, R) = ^ coll (r, R) (P S udakov(r, R)-l) + ^Lo(r, R) (4-2) 
+4>i,ho(r, R) + tppc(r, R) + A,pc(r, R) ■ 

In Eq. (|4-2p the first term represents the contribution from the Sudakov-resummed 
small-angle parton splitting; the second and third terms give the leading-order final- 
state and initial-state contributions, respectively; the last two terms come from power 
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corrections oc Qq/ Et , Qo — 2 — 3 GeV, when one integrates over the Landau pole in 
the modified leading logarithmic approximation (MLLA). This approach was shown 
to provide a very good description of the differential intra-jet energy flow at the 
Tevatronj^ 1 as measured by CDF IlP^l Thus, reliable predictions for the jet sub- 
structure in p+p reactions at RHIC and the LHC can be obtained and used as a 
baseline to study the distortion of jet shapes in more complex systems, such as p+A 
and A+A. 

Inclusive jet cross sections and jet shapes in nuclear collisions are closely re- 
latedP'S 



A °^r) -da^(R)/dy^E T J e J e 2-, (I _ (I _ Ug) . e) 3 ' 1 ->' 



X dyd?E' 



It should be noted that vacuum and medium-induced parton showers become more 
collimated with increasing E' T and the mean relative jet width, 



<M 77 S^ls) . (4-4) 



is reducedP'^ Consequently, the striking suppression pattern for jets, see for ex- 
ample the right panel of Figure [21 can be accompanied by a very modest growth in 
the observed (r/R). We show in Table Hm the relative widths in the vacuum, for a 
hypothetical case of complete parton energy loss (P 9iS (e) = 6(1 — e)) that falls inside 
of R (f q: g = 1), and for our realistic simulation of Et = 20, 50, 100, and 200 GeV 
jets of R = 0.4 in central Pb+Pb collisions at the LHC. These numerical results 
show that there is very little < 10% increase in the magnitude of this observable. 
For Au+Au and Cu+Cu collisions at RHIC, we fond that, on average, jet broadening 
is even smaller, < 5%. Therefore, a rough 1-parameter characterization of energy 
flow in jets will not resolve the effect of the QGP medium. 

Lastly, we point out where the anticipated jet broadening effects will be observed 
in the differential shape by studying the ratio iptot.( r / R)/' l Pva,c.( r / R) i n Figure HI 
The left panel shows a simulation of the differential shape of a 100 GeV jet at the 
LHC. The right panel presents the ratio discussed above. We recall that the small 
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Table III. Summary of mean relative jet radii (r/R) in the vacuum, with complete energy loss, 
and in the QGP medium. Shown are results for cone radii R = 0.4 and transverse energies 
E T = 20, 50, 100, 200 GeV at = 5.5 TeV central Pb+Pb collisions at the LHC. 




Fig. 4. Simulation of a jet shape in central Pb+Pb collisions at the LHC (left panel). The ratios 
of the medium-modified jet shape in heavy-ion collisions to the jet shape in the vacuum for jet 
energies E T = 20, 50, 100, 200 GeV for R = 0.4 (left panel). 



r/R < 0.25 region of the intra-jet energy flow in p+p collisions in our calculation has 
uncertainties associated with the normalization of the jet shape. In the moderate 
and large r/R > 0.25 region our theoretical model gives excellent descriptions of the 
Fermilab Run II (CDF II) data. The QGP effects are manifested as a suppression 
near the core and enhancement near the periphery of the jet. In the tail of the energy 
flow distribution for a cone radius R = 0.4 this ratio can reach values ~ 1.75. 



§5. Conclusions 



In summary, we discussed the discriminating power of jet observables with re- 
spect to models of parton propagation and energy loss in the QGP. More stringent 
constraints on theory, such as the ones that upcoming data on jet production and 
modification in heavy-ion reactions will provide, are necessary for precision jet to- 
mography and accurate determination of the properties of strongly-interacting nu- 
clear matter. To this end, we presented perturbative QCD calculations of inclusive 
jet cross sections/shapes and Z°/7*-tagged jets to 0[p? s ), O{a 2 s of ed ') and 0{GFCt 2 s ), 
0(GFU s (x™ ed '), respectively. To demonstrate the sensitivity of these observables to 
the characteristics of the vacuum and medium- induced parton showers, we examined 
their pronounced dependence on the cone radius R at next-to-leading order. Exper- 
imental jet physics results in heavy- ion collisions are still very preliminary but agree 
qualitatively with the theoretical expectations. Future studies in this direction will 
focus on the inclusion of non-perturbative hadronization effects and extending the 
NLO calculations of jets at RHIC and at the LHC. In their entirety, these theoret- 
ical advances will provide first-principles insights into the many-body QCD parton 
dynamics at ultra-relativistic energies. 
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